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ences  between  them. 
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PART  I 

SCATTERING  CROSS  SECTIONS  FOR  RANDOM  ROUGH  SURFACES 
—FULL  WAVE  ANALYSIS 

Ezekiel  Bahar 

1  ,  Electrical  Engineering  Department 
University  of  Nebraska-Lincoln 
Lincoln,  Nebraska  68588 

Abstract 

The  full  wave  approach  developed  earlier  to  evaluate  the  radiation  fields 
scattered  by  deterministic  two  dimensionally  rough  surfaces  is  used  here  to 
determine  the  scattering  cross  sections  for  random  rough  surfaces.  The  medium 
below  the  irregular  boundary  is  characterized  by  complex  permittivity  and 
permeability.  For  slightly  rough  surfaces,  the  full  wave  solutions  for  the 
incoherent  scattered  fields  are  shown  to  be  in  agreement  with  the  perturbation 
solution.  However,  when  the  major  contributions  to  the  scattered  fields  come 
from  the  region  of  the  rough  surface  around  the  stationary  phase  (specular) 
points,  the  full  wave  solutions  are  in  agreement  with  the  physical  optics  solutions. 
Thus,  the  full  wave  solutions  which  reduce  to  the  perturbation,  the  physical  optics 
and  the  geometrical  optics  approximations  In  special  cases,  precisely  determine 
the  limitations  of  these  approximations  and  reconcile  the  differences  between  them. 

The  full  wave  solutions  satisfy  duality,  reciprocity  and  realizability 

relations  in  electromagnetic  theory  and  they  are  invariable  under  coordinate 
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The  full  wave  solution  for  the  scattered  radiation  field  by  deterministic 
two  dimensionally  rough  surfaces  (Bahar  1980)  is  applied,  In  this  paper,  to  problems 
of  scattering  and  depolarization  by  random  rough  surfaces.  The  results  ^f  this 
analysis  are  compared  with  the  solutions  derived  from  two  general  approaches  to  random 
rough  surface  scattering  problems:  the  perturbation  technique  and  the  Kirchof f-Physical 
Optics  approximation  (Ishimaru  1978).  The  perturbation  technique,  which  applies  to 
surfaces  that  are  slightly  rough  was  used  by  Rayleigh  (Strutt  1896)  and  extended  by 
Rice  (1951),  Barrick  and  Peake  (1968),  Barrick  (1970,  1971),  Wright  (1966)  and  Rosich 
and  Wait  (1977),  Valenzuela  (1978),  The  Kirchoff-Physical  Optics  approximation  technique 
was  applied  to  surfaces  with  radii  of  curvature  that  are  much  larger  than  the  wavelength 
of  the  electromagnetic  excitation  (Beckmann  and  Spizzichino,  1963,  Beckmann  1968, 

Ament  1953). 

The  principal  elements  of  the  full  wave  approach  are  (Bahar  1980):  (a)  Complete 
expansion  of  the  fields  into  vertically  and  horizontally  polarized  waves.  The 
complete  spectrum  of  the  waves  consists  of  the  radiation  fields  (considered 
here  in  detail)  and  the  surface  and  the  lateral  wave  terms  (Bahar  1973a, b). 

(b)  Imposition  of  exact  boundary  conditions  at  the  irregular  interface  between 
two  media  y  >  h(x,z)  and  y  <  h(x,z)  characterized  by  complex  electromagnetic 
parameters  e  and  u  for  exp(iwt)  time  harmonic  excitations.  Thus  approximate 
impedance  boundary  conditions  are  not  used  in  this  work.  (c)  Use  of  Green's 
theorems  tc  avoid  term  by  term  differentiation  of  the  complete  expansions. 

(d)  Conversion  of  Maxwell's  equations  into  rigorous  sets  of  coupled  first- 
order  differential  equations  (generalized  telegraphist's  equations)  for  the 
wave  amplitudes.  (e)  Use  of  a  variable  coordinate  system  that  conforms  with  the 
local  features  of  the  rough  surface.  Thus  there  ate  no  restrictions  on  the  height 
or  slope  of  the  rough  surface,  and  both  upward  and  downward  scattering  of  the 
incident  fields  are  accounted  for  in  the  analysis.  The  effects  of  shadowing  can 
also  be  included  In  the  full  wave  analysis.  The  full  wave  solutions  are  shown  to 


satisfy  duality,  reciprocity  and  realizability  relationships  in  electromagnetic 
theory  (Bahar  1980) . 

For  the  convenience  of  the  reader,  the  full  wave  solutions  for  the  radiation 
fields  scattered  by  deterministic  two-dimensionally  rough  surfaces,  are  summarized 
in  Section  2  since  they  constitute  the  starting  point  of  the  present  analysis. 

Of  particular  interest  in  this  work,  is  the  bistatic  scattering  cross  section  per 
unit  area  for  rough  surfaces  (See  Appendix  A)  (Barrick  1970,  lshimaru  1978). 

The  full  wave  solutions  for  the  scattering  cross  sections  are  developed  first  for 
slightly  rough  surfaces  in  Section  3.  In  this  special  case  the  full  wave 
solutions  are  shown  to  be  in  complete  agreement  with  the  scattering  cross  sections 
for  the  incoherent  (diffuse)  fields  (Barrick  and  Peake  1968). 

In  Section  A  of  this  paper,  the  full  wave  approach  is  applied  to  rough  surfaces 
with  normal  height  distributions.  In  this  section  it  is  shown  that  if  a  high 
frequency  stationary  phase  approximation  is  made  to  the  full  wave  solutions,  the 
expression  for  the  scattering  cross  section  reduces  to  the  Kirchoff — Physical  Optics 
solutions,  (Beckmann  and  Spizzichino  1963,  lshimaru  1978),  with  the  exception 
that  here,  consistent  with  reciprocity,  the  Fresnel  reflection  coefficients  are 
evaluated  at  the  stationary  phase  points,  rather  than  at  the  angle  of  incidence 
with  respect  to  the  reference  plane. 

Except  for  scattering  in  the  specular  direction  with  respect  to  the 

reference  plane  (0^  “  Q^),  the  perturbation  and  the  physical  optics  solutions  are 
o  o 

not  in  agreement  even  for  slightly  rough  surfaces.  Since  both  the  perturbation 
and  physical  optics  solutions  are  derived  here  (as  special  cases)  from  the  full 
wave  solution,  the  limitations  of  each  of  the  two  special  approaches  are 
examined  and  the  differences  between  these  two  solutions  are  reconciled.  Thus 
it  is  shown  that  the  physical  optics  solutions  cannot  be  used  when  the  major 
contributions  to  the  scattered  fields  do  not  come  from  regions  of  the  rough  surface 
around  the  stationary  phase  points.  As  a  result,  the  perturbation  solution 
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(and  not  the  Physical  Optics,  solution)  should  be  used  for  slightly  rough  surfaces, 
even  when  the  radii  of  curvature  for  the  rough  surfaces  are  much  larger  than  the 
wavelength. 

For  very  rough  random  surfaces,  it  is  shown  that  in  agreement  with  the  Physical 
Optics  solution,  the  scattering  cross  section  is  proportional  to  the  probability 
density  function  for  the  rough  surface  slopes.  However,  even  at  high  frequencies 
these  solutions  cannot  be  used  when  the  incident  or  scatter  angles  are  much  larger 
than  the  mean  value  of  the  rough  surface  slope  (see  Section  A). 

The  full  wave  approach  is  not  limited  to  slightly  rough  surfaces  or  to  the 
special  cases  when  the  stationary  phase  approximations  are  valid.  However,  in  these 
two  special  cases  the  full  wave  solutions  simplify  significantly  since  they  do  not 
depend  explicitly  on  the  slopes  of  the  rough  surface.  In  Section  5  the  relation¬ 
ships  between  Physical  Optics,  Geometrical  Optics,  Perturbation  and  the  Full  Wave 
solution  are  summarized. 

It  is  interesting  to  note  that  in  order  to  obtain  the  perturbation  solution 
(as  a  limiting  case  of  the  full-wave  solution),  it  is  assumed  here  that  the  slope 
of  the  rough  surface  is  small  but  no  restrictions  are  made  on  the  height  of  the 
rough  surface.  Thus,  the  perturbatioi  solution  derived  here  contains  both  the 
Incoherent  (diffuse)  scattered  fields  (Rice  1951,  Barrick  and  Peake  1968,  Barrick 
1970,  1971,  Wright  1966,  Rosich  and  Wait  1977,  Valenzuela  1978)  as  well  as  the 
coherent  scattered  fields  (see  Section  3  and  Appendix  A) . 
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2.  Formulation  of  the  Problem 


The  starting  point  for  the  present  analysis  of  scattering  and 
depolarization  by  random  rough  surfaces  is  the  full  wave  solution  for  the 
scattered  radiation  field  by  deterministic,  two  dimensionally  rough 
surfaces  (See  Fig.  1) : 


y  -  h(x,z)  =  f(x,y,z)  =  0  (1) 

The  incident  and  scattered  radiation  fields  are  decomposed  into  a  complete 
spectrum  of  vertically  and  horizontally  polarized  components  with  respect 


to  the  reference  plane  normal  to  a^.  Denoting  the  incident  and  scattered 


fields  by  the  superscripts  i  and  f  respectively  and  the  vertically  and  horizontally 
polarized  components  of  the  fields  by  the  superscripts  V  and  H  respectively, 
the  electric  and  magnetic  fields  E  and  II  can  be  expressed  in  matrix 
notation  as  follows: 


,4 


where  n  =  (u  /e  )  is  the  intrinsic  impedance  for  free  space.  The  full 
o  o  o 


,f  . 


wave  solution  for  G  is  (Bahar  1980) 
„f 


ClnTfFT1  exp[ik  (nf-n1)'?  JU (7  )dA*n  G1 
o  r  o  s  s 


2  Go  C(nf,ni)G1 


(3) 


An  exp(iwt)  time  dependence  is  assumed  and  the  constant  G  is  given  by 


G  *  -ik  exp‘L-ik  r^j/2irr* 
o  o  o 


(4) 


... 


where  k  *  l>(u  /c  )  a  is  the  free  space  wave  number  and 
o  o  o 


rf  f  ~f  f/  •  nf  J  -  Qf  -  ,  Qf  .  ,f  -  , 

r  *  r  n  *  r  (sinG  cosd)  a  +  cos0  a  +  sino  sind)  a  ) 
o  x  o  y  o  z 


(5) 


is  the  position  vector  to  the  observation  point.  The  position  vector 


>  \ 


EVi 

HVi 

HVf' 

G1  - 

=  no 

,  g£  - 

■  n 

.  (2) 

EHi 

o 

HHi 

E»f 

0 

HHfJ 
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to  the  source  is 

r1  =  -r1  n1  =  -r1[sin01  cosdi1  a  -  cosG1  a  +  sinG1  sinrf)1  a  )  !  > 

o  r  x  o  y  o  v  z' 

in  which  n1  is  the  direction  of  propagation  of  the  incident  waves. 

The  position  vector  to  the  rough  surface  is 

r  =xa+h(x,z)a  +za  =r-f(x,y,z)a  , 
s  x  y  z  y 

in  which  f(x,y,z)  is  given  by  (1)-  The  elementary  area  of  the  rough 

surface  is 

dA  =  n  dx  dy / (n-a  )  . 

y 

in  which  n  is  the  unit  vector  normal  to  the  rough  surface: 

n  =  Vf /  I Vf  I  =  [-h  i  +  a  -  h  a  J/(h2  +  1  +  h2)*5 
11  xx  y  z  z  x  z 


-  sinv  cosfia  +  cosy  a  +  sinv  sin6  a 
x  y  1  z 


(6) 


(7) 


(8) 


(9) 


and 


h^  =  3h/3x  ,  hz  =  3h/3z 


(10) 


A  local  Cartesian  coordinate  system  with  coordinate  surfaces  normal  to 
the  unit  vectors  and  is  employed  to  derive  the  full  wave  solution: 

n  =  n  x (a  x  n)/|a  x  n|,  n„  =  n  and  n  =  (eu  x  n)/|a  x  n|  (11) 

X  X  X  Z  J  X  1 

The  vertically  and  horizontally  polarized  components  of  the  incident  and 
scattered  electric  and  magnetic  fields  with  respect  to  the  local  plane, 
normal  to  the  unit  vector  u ,  are  denoted  by  the  subscript  n. 

They  are  related  to  the  components  with  respect  to  the  reference  plane 


through  the  transformations 


fEVi' 


.in  _  | 


l, Hi 

iLn 


T1  G1 


r  <  Si1 

i|)  \jj 

-si  c!i 

i pj 


i^  r_Vi 

E 


.Hi 


(12) 
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J.  ,.fn 
=  T  G 


fcf  -sf 1 iEv£i 


in  which  C  and  s|  are  the  cosine  and  sine  of  the  angle  between  the  local 
V  v 

plane  of  incidence  and  reference  plane  of  incidence  normal  to  the  unit  vectors 
a^  and  a  » respectively .  Thus  they  can  be  expressed  in  terms  of  the 
scalar  product  and  the  scalar  triple  product: 

=  C°S^  =  SHi*SHi  *  =  Sin^  =  [SHi5Si"1]  ’  (1 


/  i  \  / 1  —  i  i  “  n  -i  ■*" .  / 1  “  x  ~  i 

aHi  =  (n  x  ay)/ln  x  a  |  ,  aHi  =  (nxn)/|nxn| 

Similarly  C  and  S,  are  the  cosine  and  sine  of  the  angle  between  the  local 
plane  of  scatter  and  the  reference  plane  of  scatter  normal  to  unit  vectors 
and  a^  .respectively .  Thus 


_f  .  f  -  -n  „f  .  ,f  -n  -fn 

S  =  C°SV  =  aHf‘aHf  ’  =  Slnlp  =  LaHf  aHf  n  ] 


where 


aHf  =  (°fx  ay)/lnfx  ayi  ’  aHf  =  (°fx  n)/|nfx  n| 

The  elements  FP<^(P,Q=V,H)  of  the  2x2  local  scattering  matrix  F(nfn,nin) 
in  (3)  are  (Bahar  1980): 


CinFVV 

o 


2C^?\(y  C^fncos((|)f-  A-sV")  (1-1/e  )+(l-u  )cos(<t!f^  f)j 

\cl"><C+\c^<C+co£n) 

2C^  (c,C?tfncos(*fn-r)-Sl®  (l-l/hj+(l-e_)cos(0f%in)  j 


..in^Hii  o  o  r  1  1 

C  F  «*  - - - - 


(coin</nr)(cf;^/nr)(c^cofn) 


CinF1,V 

O 


CinFVH 

o 


-sin  (f.b‘-<T)2C^nri  (l-l/tr)cJ-(l-l/pr)C*'] 
<+nrc;b)(cof^r)(c>c^) 
sin(^-^;2C*\^n  [(1-1/p  )C^(l-l/e  )<??. 

O  O  L  L  J.  r  a. 

/„in_in,  ,  fn,  fh  in  ,  ,,fn. 

(Co  *lci^rir )  (CQ  +  Co  ) 


(18d) 
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in  which  Che  dimensionless  quantities  ,  r]  ,  t  and  are  the  refractive 
index,  relative  intrinsic  impedance,  relative  permittivity  and  relative 
permeability  ,respec  tively , 

nr  =  .  nr  =  -:i/r,o  =  (VV 

1  o 

c  =  e./c  ,  V  =  p./p  .  a 

r  1  o  r  t  o 

The  media  for  y  >  h(x,z)  and  y  <  h(x,z)  are  denoted  by  the  subcripts  o  and  1 
respectively.  The  permittivity  and  permeability  for  medium  1  can  in  general, 
be  complex  to  account  for  dissipation.  The  cosines  and  sines  of  the  angles  of 

incidence  and  scatter  (with  respect  to  the  local  coordinate  system)  0in  and 

o 

0^n  in  free  space,  v  >  h(x,z),  are  given  by: 

-in  m  -i  -f  fn  -fn  -f  -  ,, 

C  =  coso  =  -n  *n  ,  C  =  cos0  =  n  *n  >  (2 

o  o  o  o 


in  m 

-l  -f  ,,fn 

=  cosu 

-n  *n  ,  C 

o  o 

0 

,in  i-i 

„fn  ; — f 

'  =  In  x  n 

,  S  =  n 

o  1 

o  1 

The  sines  of  the  corresponding  angles  in  medium  1,  y  <  h(x,z),  are  given  by 
Snell ' s  law  • 

cin  .  .in  _in,  fn  .  „fn  „fn. 

S  =  sin0  =  S  /n  ,  S,  =  sin0,  =  S  /n 
1  lor  1  lor 

Thus 

„in  -in  r  in  2-,^  fn  fn  .cfn.2,ii 

C1  =  cosO^  =  Ll-CSj^  )  j  ,  =  cosOj^  =  Ll-(S1  )  J 

The  cosine  and  sine  of  the  angle  between  the  planes  of  incidence  and  scatter 

in  the  local  coordinate  system,  (1)  are  given  by: 

„ , fn  , in.  -n  -n 
c°s«  )  =  aHf-aHj, 

and 

,  .  fa  in  , -n  -n  - 

sin(-J)  -<f  )  =  iaHf  a}(.  nj  . 

The  shadow  function  L'(r  )  in  (3)  is 


)-  -I 


1  ,  illuminated  and  visible  region 
(o  ,  nonilluminated  or  nonvisible  region 


1-8 


The  nonilluminated  (shadow)  region  of  the  rough  surface  extends  from  the 

locus  of  points  r  =  r1.  to  the  locus  of  points  r  =  r1  that  satisfy 
s  si  s  sz 

ni*n(ri.)  =  0  and  (r*  -  ri.)*n(ri,)  =  0  •  (25) 

si  s2  si  si 

Similarly,  the  nonvisible  region  of  the  rough  surface  extends  from  the  locus 

of  points  r  =  r£,  to  the  locus  of  points  r  =  rr_  that  satisfy 
r  s  si  r  s  s2 

n£*n(r£  )  =  0  and  (r£  -  r^1)*n(r^1)  0  •  (26) 

si  s2  si  si 

The  full  wave  solution,  (3)  satisfies  reciprocity,  duality  and  realizability 
relations  in  electromagnetic  theory  (Bahar  1980).  Of  special  interest  in 
this  work  is  the  normalized  scattering  cross  section  per  unit  area  for  rough 
surfaces  (Ishimaru  1978): 


oP<)  .  4,(r£)2iEP£|2/Ayi,1|2  .  < 

in  which  A^  is  the  projection  of  the  area  of  the  rough  surface  A  on  the 
reference  plane  normal  to  a^.  Thus  for  P,Q  c  V,H 

A  lcpy.  A  jDpy  ..pUk  -?;)]»(? 


in  which  the  symbol  *  denotes  the  complex  conjugate , and  D 


(n*a  )(n'*a  ) 

y  y 
(28) 


are  the  elements  of  the  matrix  D: 

|C£ (FVVC1-FVHSi)"Sf (FHVCi-FHHS1)  Cf (FVVS*+FVHC1)-S£ (FHVs|+FHHC1)  | 
D  =  Cln  T£FTl=  Cln| 

°  °  !sf(FVVCi-FVHSi)-K:f(FHVC£-F1,HS£)  Sf  (F^+F^f  )4Cf  (F^sf+F^f ) 

[_  v  o  V  y  <J>  V  v  'P  r  <r  'P  <P  J 


Since  D^yr  )  is  a  function  of  the  angles  of  incidence  and  scatter  in  the 

local  coordinate  system,  (11),  it  depends  explicitly  only  on  the  gradient 

of  the  rougli  surface  Vf  (9).  The  exponent 

e  x  p ;  ’ k  (n£-n1)-(r  -r  '  )  J  -  exp[ iv  (x-x')+  iv  (z-z')+  iv  (h (x , z)-h (x ' , z 1 ) j  ,  (30) 
o  s  s  x  z  y 
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depends  explicitly  only  on  the  height  o£  the  rough  surface  h(x,j).  The 
vector  v  defined  in  (6.30)  is  given  by 

v  ■=  v  a  +  va  +  v  a  =  k  (n*1-  nl)  .  (31) 

xyyyzzo 

The  shadow  function  U(r^)  (24)  in  general,  depends  on  the  gradient  of  the 

rough  surface,  through  the  unit  vector  n  (9),  as  well  as  on  the  height  of 

the  rough  surface  h(x,y),  (25),  (27).  However  L^r^)  /  1  only  if  n^n  >  0 

or  n^'n  <  0  on  portions  of  the  rough  surface  independent  of  h(x,z).  Thus, 

the  shadow  function  U(r  )  is  more  sensitive  to  the  gradient  of  the  rough 

s 

surface  than  to  its  height (Saneer  1969,  Brown  1978). 

The  remainder  of  this  paper  deals  with  random  rough  surfaces  for  which 
only  the  statistics  of  the  height  h(x,z)  and  its  gradient  are  assumed  to  be 
known.  Section  3  deals  with  slightly  rough  surfaces  and  the  full  wave 
solutions  are  compared  with  the  perturbational  solution  (Rice  1951,  Barrickand  Peake 
1968,Bairick  1970,Rosich  and  Wait  1977).  In  section  4,  h(x,z)  is  assumed  to  be  normall 

distributed  and  no  restrictions  are  made  on  the  variance  of  the  rough  surface. 

These  solutions  are  compared  witii  both  the  perturbational  and  physical  optics 
solutions  (Beckmann  and  Spizzichino  1963). 


3.  Incoherent  Scat.ering  Cross  Section  Per  Unit 
Area  For  Slightly  Rough  Surfaces. 

For  the  incoherent,  diffuse  field,  the  scattering  cross  section  for  unit 
cross  sectional  area  j.s  given  by  (Ishimaru  1978)  : 

<  aPQ  >  =  4TT(r1)2  <  |Epf-<Epf>  j2  >  /Ay|EQl|2  ,  (32) 

in  which  the  syiui.ox  <  >  denotes  the  statistical  average.  For  slightly  rough  surfaces 
tan  y  <<  1  and  n  is  set  equal  to  a^  in  (3).  Thus 


T1  -  I  ,  Tf  ->  I  ,  U(r  )  -  1  ,  n  •  a  -  1 

s  y 


(33) 


in  which  I  is  the  2x2  identity  matrix.  Thus,  the  elements  of  the  matrix  D  (29) 


.me  no  longer  functions  of  position  and  may  be  extracted  from  the  integral,  (32). 
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The  incoherent  scattering  cross  section,  (32)  reduces  to 


PQ>  ‘XQ| 

<  a  H  >  =*  - 

O  TT 


■3  . 


(34) 


in  which 


upQ  =  (cinFPQ^_  _ 
o  o  n=a 


(35) 


and 


ert  «=  — 

J  A 

y 


2- 

expliv  (x-x)+  iv  (z-z ’ ) ][x, (v  ,-v  )-jX(v  )]  _dxdy  dx'dy'  .  (36) 

x  z  z  y  y  y 


The  characteristic  function  x(vy)  and  c*ie  joint  characteristic  function 
v(v  ,-v  )  are  defined  in  terms  of  the  probability  density  function  W(h)  and 

y  y 


joint  probability  density  function  W(h,h'): 


X(v  )  =  J  exp(iv  h)W(h)dh 
y  )  y 


and 


X(vyi-vy)  =  ]  exp[iv  ;(h-h' ) ]W(h,h' )dh  dh' 


(37) 


(38) 


Thus  for 


and 


slightly  rough  surfaces 


exp(iv  h)  =  1  +  iv  h  -  h^/2 

y  y  y 


(39) 


12  2  12  2 

X(v  )  -  1  -  —  v  <h  >  =  1  -  —  v  a  ,  (40) 

y  2  y  2  y  o 

in  which  0  is  the  variance  of  h  and  <h>  *  0.  Similarly, 

o 

X.,  (v  ,-v  )  -  1  -  v^  <  h^  >  +  <  hh'  >  =  1  -  v^  a  (1-C)  ,  (41) 

^2  y  y  y  y  y  o 

in  which  C  is  the  normalized  correlation  coefficient.  Assuming  that  the  rough 

surface  is  statistically  homogeneous  and  isotropic,  the  surface  height 

correlation  function  <  hh '  >  depends  only  on 

x,  -  x-x'  and  z,  -  z-z'  .  (42) 

d  d 

Thus  if  the  correlation  distance  is  auch  smaller  than  the  width  of  the  illuminated 


surface,  (36)  reduces  to 
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&  "  expLiv  x  +  iv  z.J  <  hh’>  v  dx_,  dzJ 
J  \  xd  z  d  ydd 


2  2 

it  v  W(v  ,v  )  . 

y  x  z 


(A3) 


in  which  w(vx»vz)  is  the  spectral  density  of  the  rough  surface  height 
function  (Barrick  1970,  Ishiraaru  1978).  Thus  the  incoherent  scattering 
cross  section  for  slightly  rough  surfaces  is  given  by 

2 

n-+a 


PQ  4 

<0  >  =  TtkW(v,v) 

o  o  x  z 


ClnFPQ(cosGf  +  cosQ1) 
o  o  o 


(44) 


The  expression  (44)  is  precisely  equal  to  the  perturbation  solution  for  the 
scattering  cross  section  for  slightly  rough  surfaces  (Barrick  and  Peake 
1968,  Barrick  1970).  In  the  notation  used  for  the  perturbation  solution 


[C  F  ^(cos +  cosO^)]-  -  =  t  2  cos0P  cos0i  a1  ^ 

o  o  o  n-*a  o  o  * 

y 

in  which  the  upper  and  lower  signs  are  used  for  Q  =  H  and  V,  respectively. 
For  the  perfectly  conducting  case  the  scattering  cross  section  is  given  in 
matrix  form  by 


(45) 


<C  >  =  4tt  k  W(v,v) 

o  O  X  z 


[sinQ1  sin©1  -  cos  (ij^-ij)1)  ]  [cosU1  sin($*-$*)  V 
o  o  o 


f  f  i  2 

[cosOq  sin(4>  — )  j 


[cosO^cosG^cos  (4)^ -i>1)  j 


For  backscatter  0^  =  0*,  ij)  —  4> 1  =  u  (46)  reduces  to 


0  >  =  4tt  k^  W(2k  sin0i,  0) 

0  B  o  o  o 


2  i  2 
(sin  6  +  l)z 

o 

0 


4  l 
cos  0 

o 


(47) 


4.  Incoherent  Scattering  Cross  Section  For  Rough 
Surfaces  With  Normal  Distributions 


When  the  slope  of  the  rough  surface  is  not  small,  the  incoherent 
scattering  cross  section  for  high  frequencies  ,  (32)  can  be  written  as  follows 
(see  Appendix  A)  : 


(46) 
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! 

1 


i 


k2 , 

FQ  O  »  PQ  2 

<  0  >  -  —  j  '  S  '  ^2(vy,-Vy>^'X(vy)(  >xptivxxd+  ivasxd.dxd  dzd  . 

A 

y 

in  which 

„  DPQ(?  )D?Q*(r’)U(F  )U(r')  i DP V  )U(?  ){2 

,  ccQ  a  s  s  s  I  t>  r  1 


(n*a  )(n'*a  ) 

y  y 


(n»a  ) 

y 


(48a) 


(48b) 


To  derive  (48)  it  is  assumed  chat  the  surface  height  and  gradient  are  uncorrelated  (whici. 

is  true  for  Gaussian  surfaces)  and  that  the  probability  density  function  tor  the  shadow 

function  U(r^)  depends  only  upon  the  gradient  of  the  rough  surface  (see  Section  2> 

As  in  Section  3,  the  rough  surface  is  assumed  to  be  statistically  homogeneous 

po 

and  isotropic*  thus  *  S  >  and  %  depend  only  on  x,  and  z  ,  (42). 

^  d  d 

For  normally  distributed  surfaces  the  probability  density  function  and  the 
joint  probability  density  functions  W(h)  and  U(h.h')  are  given  by: 


w  (h)  “  - T -  exp(-~ ) 

(2ir)  J  <J  2-r 


(tV) 


and 


W(h.h’) 


2 ft  .J2(l-C2)‘ 


exp  - 


,  h  -  2Chh'+  h * 

2  2 
2o Cl-C ) 


(50) 


in  which  0q  is  the  variance  and  C  is  the  correlation  coefficient.  The  corresponding 
characteristic  functions  are 


and 


X(vy)  -  exp(-v2  j2/2) 

2  2 

X2<vy.-v  )  -  exp[-vz  Op(l-C)i 


(51) 

(52) 


Thus  it  follows  that 


'•2<Vy’~V  ~  h(vy)|2  "  exp(-v2  o‘)[fcxp(v2  o2  C)  -  1  j 

and  the  sui face  height  correlation  function  is 

,  2 

•  h(x,y)  tli  (x'  ,  y ' )  >  u  C 


(53) 


(54) 


I 

I 
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Thus  for  slightly  rough  surfaces  f(v  aQ)  «  1,  n  -  a  "[  ,  (54)  reduces  Co 


X2(Vy*“V  "  !x(vy)|2  -  vy  <  hh’> 


<  SP^  >  =  [SP^ 


-  -  -  !dpq|2 

n-*a  '  o  1 
v 


pQ 

and  <  u  H  >  is  given  by  the  perturbation  solution,  (44), 

II  it  is  not  assumed  apriori  that  the  surface  is  slightly  rough,  *•  a 
can  be  expressed  approximately  as  follows  for  high  frequenc ies: 

71  »  ./ 


PQ 


In  which 


ffm  2n 


Jo(Vxzp)exp("Vv  ao^exp(vv  °o  c)  "  1-:'pdp 

oo  J  J 


where  J  is  the  Bessel  Function  of  the  first  kind,  and 
o 


,  z  ,  2  ,  .  i  .  z. 

vxz  =  (vx  +  V  and  p  =  (xd  +  V 

The  correlation  coefficient  C  is  only  a  function  of  p. 

Assuming  a  correlation  coefficient  of  the  form 

C(p)  =  exp(-p2/£.2)  ,  (58) 

where  £  is  the  correlation  distance,  it  can  be  shown  that  (Beckmann  and 
Spizzichino  19b3): 

[u(v  a  £)2exp(-v2  o2  -  v2  /4)  =  tt2  v2  w(v  ,v  ),(v  a  )2  << 
y  o  r  y  o  xz  y  x’  z  y  o' 

J  J  # 2  ,2  <Vy  °o)2m  .2.2... 

^  =  y£  exp (-V  oq)2:  — -  exp(-vxz  £  /4),  v  0Q  =  1. 

|  m=l 

II  £2  r  ,  ,  ,2-|  ,  ,2  ..  , 

— z - t~  expL-  (v  £/2v  o  )  j  ,  (v  O  )  1  • 

2  2  xz  vo  V  o 

.  v  a 
(y  o 

PO  PO 

In  (56)  S  y  is  .:ed  by  its  value  at  the  stationary  phase  points  of  the 

integrand  in  (3).  Thus  (Bahar  1980): 

SIQ  -  £co«0ol  F3(Ot>;,0f-,fri)j2(RPo(o“)!26pQ  . 


2^ 


(55a) 

(55b) 


(56a) 

(56b) 

(57) 


1, 

(59) 


(60) 
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PQ 

In  view  of  the  Kronecker  delta,  6^,  Sc  vanishes  for  P#Q. 

g 

In  (60)  0^  is  the  local  angle  of  incidence  and  scatter  at  the  specular 

points.  Thus,  , 

cos  0S  ■*  (1-n  •n1)/2 
o 

“tit  cosG*  cosG^  -  sinO'*'  sinO^  cos(<p^-tp^)  j/2  ,  (61) 

o  o  o  o 

P  s 

The  Fresnel  reflection  coefficients  R,  (0  )(P=V,H)  are  evaluated  at  the 

lo  o 

^  £  g  J 

stationary  phase  points  where  0  ►  0  >  (j  and  not  at  the  angle  of  incidence  0  as  done 

'  ooo  o 

inBeckmann  and  Spizzichino , 1963 .  ( Ishimaru ,  197?).  The  coefficient  is  precisely 

the  expression  derived  for  the  Physical  Optics  solution  (Beckmann  and  Spizzichino 
1963) 

iff!  ^  ou=  On 

(62) 


f  (eSe^V)  = 

J  o  o 


2  cos2  0s 
o 


cosG'*'  (cosG^  +  cosG"”") 
o  o  o 


PQ 

Thus  the  high  frequency  approximation  for  <  a  >  is  given 


by 


PQ  ,  o 
a  >  ~  — 

a)  n 


2,s 


2  cos  6 


cosG^  +  cosG'*’ 
o  o 


i<X>i2 


PQ 


(63) 


P  PP 

For  the  perfectly  conducting  case  jRp0l  =  1  and  <  Om  >  is  independent  of 
polarization.  Thus  the  full  wave  solution  for  the  incoherent  scattering  cross 
section  reduces  to  previously  derived  Physical  Optics  results  when  the  following 
conditions  are  satisfied. 

(a)  For  high  frequencies,  the  principal  contributions  to  the  incoherent 
scattered  field  come  from  the  stationary  phase  points  of  the  integrand  in  (3), 
(flaliar  1980).  In  this  case  n  can  be  replaced  by  its  value  ng  at  the  specularly 
oriented  portions  of  the  rough  surface  : 


n^  «  (n^-n*)/ |  n^-n^  | 

(b)  The  vectors  n^  n^  and  a  are  coplanar.  Thus 

“f  -  -  .  ;  “I  -  -  •  .  — i  —  — f  |  . 

n  a  n  j  =  n  a  n  j  =  in  a  nj  =  0 
>'  s  y  sJ  y 


and  at  the  stationary  phase  points 


.  f  |  i 

ip  -  <P 


(64) 


(65) 


(66) 
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In  this  case  the  local  and  reference  planes  of  incidence  and  scatter  are 
identical: 


T1 


Tf  =  I 


(57) 


where  I  is  the  identity  matrix 

(c)  The  phase  of  the  integrand  in  (1)  varies  rapidly  for  values  of 

r^_  in  the  shadow  region,  thus  the  value  of  the  integral  does  not  change 

significantly  on  replacing  U(rg)  by  unity  everywhere  on  the  surface. 

For  slightly  rough  surfaces  (63)  reduces  to 

PP  4 ,  p  c  ,2  a  e 

<  0  >  «  4 irk  R,  (6  )  cos  8  W(v  ,v  ) 

00  o  lo  o  O  X  z 


(68) 


The  above  expression  which  satisfies  reciprocity,  is  in  agreement  with  the 
Physical  Optics  solution  given  by  Ishimaru  (1978)  with  the  exception  that  here 
the  reflection  coefficient  is  evaluated  at  the  specular  angle^  as  required  by  the 
stationary  phase  condition)  and  not  at  the  angle  of  incidence.  For  backscatter 
(n  =  -nx  =  ngl  cos0o  =  1)  (68)  reduces  to 

<  °lP  >B°  4TrkolRlo  <0)|2w(2  ko  sin0o>  0)  •  (69) 

Note  that  for  slightly  rough  surfaces  the  perturbation  solution  (44)  is  not 

in  agreement  with  the  Physical  Optics  solution,  (68)  except  at  normal  incidence. 

This  is  because  even  at  high  frequencies  condition  (a)  is  not  satisfied  for 

slightly  rough  surfaces  except  near  normal  incidence.  For  large  angles  of 

incidence  .there  are  no  stationary  phase  (specular)  points  on  the  slightly  rough 

surface,  thus  even  for  high  frequencies,  the  Physical  Optics  solution  cannot, in 

general,  be  used  for  slightly  rough  surfaces  when  is  smaller  than  the  mean 

value  of  the  slope  f  the  rough  surface  3Q> where 


tang  =  2o  / i 
o  o 


(70) 


For  very  rough  surf aces ,[  (v^  U^)  >>  lj,the  scattering  cross  section  is  given  by: 


<o„  >  =  —  — t~ — rr  tx^(v„l/2 

o  coso  +  cosO 


!v  u  ) 
V  O 


-  if  aecHY  P(h  ,h  )|RR  (0s)  1 2  ,  ( 

s  xsp  zsp  'loo' 

in  which  p(h  ,h  )  is  the  joint  probability  density  function  for  the  slopes, 

AoD  Zbp 


(^),  at  the  stationary  phase  points  : 


p(h  ,h  )  =  - •=■ 

xsp  zsp  2  a2 

s 


h2  +h2  ] 

xsp  ZSP! 


i  2 

|  Lan  y 


-exp- 


sec^y  *  1  +  tan^y  *  l/(a  *n  =  [  2cos0S/ (cosO^  +  cos0*)'l^  a ^  *  2 .  (71 
s  s  y  s  o  o  o  s  o 


s  i  f 

Thus  for  backscatter  (Q  =  0,  <p  =  0,  <j)  =  n),  (71)  reduces 


iR^(°)|2 


tan^G* 

o 


Vi  l<o(0)|2p(ta„0;,0)  , 


<“">»  *  --o ---c-  ”■>  -  7-57  -  -Vi  l»L«»i2’““«>  •  <”> 

cos  0  tan  3  tan  3  cos  0 

p  o  0  (_  oj  o 

in  which  R,  (0)  is  the  Fresnel  reflection  coefficient  for  normal  incidence, 
lo 

Thus  the  full  wave  solution}  (43)  is  in  agreement  with  the  earlier  high  frequency 
results  (Uarrick  and  Peake  1968).  However,  for  very  rough  surfaces  the  Physical 
Optics  solution,  (72)  cannot  be  used  near  grazing  angles,  even  for  very  high 
frequencies,  since  , in  these  cases  assumptions  (a)  and  (c )  are  not  satisfied , and 
shadowing  effects  become  important. 

For  the  Physical  Optics, or  the  perturbation  approximations, derived  in 
sections  3  and  4,  it  is  not  necessary  to  know  the  complete  statistics  of  the  gradient 
Vf  of  the  rough  surface.  However,  when  these  approximations  are  not  valid,  it  is 
necessary  to  determine  the  probability  density  function  of  the  gradient  in 
order  to  evaluate  the  scattering  cross  sections  <o^>  (48). 


5.  Relationships  Between  Physical  Optics,  Geometrical  Optics, 
Perturbation  Theory  and  the  Full  Wave  Approach 

For  very  high  frequencies  (3)  can  be  written  as 


G  =  G  D(r  )U(r  )  exp[iv*r  j -  >  (73) 

o  s  s  -  s  -  - 

L  3n  n  *a 

a  s  s  y 

where  n  is  replaced  by  ng >its  value  at  the  stationary  phase  points  (60), 

v  is  given  by  (31)  and  G  is  defined  by  (4)*  Equation  (73)  is  the  physical 

optics  approximation  of  (3).  To  obtain  the  corresponding  geometrical  optics  approxir 

tion,  start  by  expanding  v-rs  about  its  value  at  a  stationary  phase  point  (where 

r  =  r  )  •  Thus 
s  so 

v-r  “vx+v  z+v  [h+h  (x-x  )  +  h  (z-z  ) J 
sx  z  yoxoo  zoo 

+  — ^r[h  (x-x  )2+  2h  (x-x  )(z~z  )  +  h  (z-z  )2I 

2  xxo  o  xzo  o  o'  zz'  o  J 


v  p 

v-r  +  — ^  h  (x-x  )  +  2h  (x-x  ) (z-z  )  +  h  (z-z  ) 
so  2  xxo  o  xzo  o'  o  zz  o' 


in  which 


h  -  (3h/3x)~  ,  h  =  (3  h/3x  )-  ,  h  =  (3  h/Jxdz)-  , 

xo  r  xxo  r  xzo  r 

so  so  so 

h  =  (3h/'Jz)-  and  h  =  (B'di/Sz^)-  , 

zo  r  zzo  r 

so  so 


v+vh  =0,v+vh  =0  .  (76 

x  y  xo  z  y  zo 

Using  a  principal  axis  coordinate  transformation  about  the  stationary  phase  point, 

r  »  the  geometrical  opt  iCs  contribution  from  the  neighborhood  of  this  point  can  be 
so 

expressed  as  (Barrie  ’970) 

i - i  _  n  it  fl  \  i t  / \  _ 1  J ~  —  1 


-  G  D(r  )U 

8° 

o 

so 

~iv 

* 

exp 

_ Z£ 

2 
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G  D(r  )U(r  )expj"ivr 
o  so  so  L 


soJ  v„„  xp  2p 


yp 


=  (Z2i)Go  expTiv 

O 


•r  1 
soJ 


■Ip  2p 


<  <o“> 

lo  o 


0 


r“  (6s) 

lo  o 


U(r  )  , 
so 


(77) 


The  integrals  in  (77) 
(Abramowi tz  and  Stegun  1964)  and  r 


are  identified  with  the  Fresnel  Integrals 
-1  -1 


lp  h  *  2p 
xxp 


are  the  principal  radii  of 


zzp 


curvature  (convex  side  facing  observer)at  the  stationary  phase  point, (v  ^  =  2k^cos6o) 

p  s  s 

and  R,  (0  )  is  the  Fresnel  reflection  coefficient  for  the  specular  angle  0  .  (65)  . 
lo  o  o' 

Provided  that  this  stationary  phase  point  is  visible  and  illuminated  [U(r  )  =  1] 
its  contribution  to  the  scattering  cross  section  (28)  is  given  by  : 


PP 

J 

8° 


it  r ,  r 
IP  2p 


R^  (0s) 

lo '  o 


(78) 


To  determine  the  average  scattering  cross  section  per  unit  area  assume  that 
the  principal  radii  of  curvature  (r^.r^at  the  stationary  phase  (specular)  points 
and  the  location  of  these  points  are  independent  random  variables.  Furthermore 
for  very  rough  surfaces  the  phase  v*r  can  be  assumed  to  be  uniformly  distributed 
from  -it  to  it,  (Beckmann  and  Spizzicnino  1963).  Thus  it  follows  from  (77)  that 

l  2 


PP  P  s  -f  -i  , - 

<0  >  =*  TT<r  r  N>  R  (0  )  P  (n  ,n  n  ) 

go  lp  2p  1  lo  o  1  2  *  1  s 


(79) 


in  which  <r^pr2p^>  Is  Lhe  avera8e  01  the  product  of  the  principal  radii  of  curvature 

and  N  the  number  of  specular  points  per  unit  area.  The  probability  that  a  specular 

point  (with  slope  n  *>  n  )  will  be  both  illuminated  and  visible  is  given  by  P„(n^,n*jn  ), 

s  z  s 

The  expression  for  P9  has  been  given  by  Sancer  (1969)  for  rough  surfaces  with 
normal,  distributions.  The  physical  optics  result  derived  by  Kodis  (1966)  for 
perfectly  conducting  surfaces  is 
PP 

■  a  •  =  tt  sr ,  r„  >  <N> 

lp  2p 


(80) 
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Barrick  (1968)  determined  the  relationships  between  the  surface  statistics 

and  <r,  r  N>  and  has  shown  that  (80)  reduces  to 
lp  2p 

PP  4 

<0  >  =  *sec  YsP(hxsp,hzsp)  ,  (81) 

which  is  in  agreement  with  (71).  In  the  derivation  of  (81),  however,  it  is 

net  necessary  to  assume  that  the  rough  surface  is  normally  distributed. 

Furthermore.it  is  incorrect  to  assume  as  implied  by  (SO)  that  r.  r„  and  N 

lp  2p 

are  statistically  independent. 

It  has  been  shown  that  in  order  to  obtain  the  physical  optics  solution  (that  has 

appeared  previously)  from  the  full  wave  solution,  the  coefficient  of  exp[iv-r  j  in  the 

integrand  of  (3)is  replaced  by  Its  value  at  the  stationary  phase  points  (n=n  ) (Section  4). 

s 

Furthermore,  to  obtain  the  geometrical  optics  solutions,  the  full  wave  integral  of  (3) 
is  evaluated  using  the  stationary  phase  (or  steepest  descent)  method.  It  has  also 
been  shown  (Section  3)  that  for  slightly  rough  surfaces,  the  full  wave  solution 
reduces  to  the  perturbation  solution  upon  replacing  the  unit  vector  normal  to  the 
rough  surface,  n,  by  the  unit  vector  normal  to  the  reference  surface,  a  ,  in 
the  coefficient  of  exp^iv’r^],  (3). 

6.  Concluding  Remarks 

The  full  wave  approach  has  been  applied  to  problems  of  scattering  and 
dcpolarizatlt  a  of  radio  waves  by  random  rough  surfaces.  The  general,  full  wave, 
expression  for  the  scattering  cross  sections  per  unit  area  is  given  in 
Appendix  A  for  the  incoherent  and  coherent  fields.  In  Section  3, which  deals  with 
slightly  rough  surfaces,  the  full  wave  solutions  are  shown  to  reduce  to  the 
perturbation  solution.  For  high  frequencies  it  is  shown  that  when  the  stationary 


phase  approxim.it  ion  is  valid,  the  full  wave  solutions  reduce  to  the  physical 
optics — Kirchoff  solution.  Thus  the  two  general  approaches  applied  to  scattering 
by  random  rough  surfaces,  perturbation  and  physical  optics,  are  derived  here  as 
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special  cases  of  the  full  wave  solution.  The  full  wave  approach  precisely 
determines  the  limitations  of  the  earlier  approaches  as  well  as  reconciles  the 
differences  between  them.  Thus,  physical  optics,  geometric  optics  and 
perturbation  theory  areallspecial  cases  of  the  full  wave  approach  (Sec.  5). 

The  full  wave  solution  which  is  invariant  under  coordinate  transformations  also 
satisfies  reciprocity,  duality  and  realizability  in  electromagnetic  theory 
(Bahar  1980) . 
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Appendix  A 


Substituting  the  full  wave  solution  (3)  into  (32)  and  assuming  that 
the  rough  surface  height  and  gradient  are  uncorrelated  and  that  the  probability 
density  function  for  the  shadow  function  U(r  )  depends  only  upon  the  gradient  of 
the  rough  surface,  the  expression  for  the  incoherent  scattering  cross  section  reduces  to 

?  OO 


:aP(3>  - 


7T 


<SPQ>  X2(vy,-vy)- 


A  X(v  , 


exp[  iv  x  +  iv  z  ,  ]dx  .  dz, 
x  d  '/.d  d  d 


(A .  1 ) 


For  simplicity  the  rough  surface  is  also  assumed  to  be  statistically  homogeneous 
and  isotropic,  thus  the  surface  height  correlation  function  <h(x,z)h(x' , z ' ) > 
depends  only  upon  the  distances 

r  =  (x-x')a  +  (z-z’)a  =  x  a  +  z  a  (A. 2) 

d  x  z  d  x  d  z 

In  (A.l)  x  anc*  X2  are  the  characteristic  function  (37)  and  the  joint  characteristic 
function  (38).  Furthermore, 

SPQ  =  DPQ(?)DPQ*(?')U(?)U(?’)/(n-ay)(H’-ay)  (A-3) 

PQ 

in  which  D  is  given  by  (29)  and  U  is  the  shadow  function  (24).  The 

symbol  <  >  in  (48)  denotes  the  statistical  average,  thus  in  order  to  evaluate 

PQ  PQ 

<S  x>  for  the  general  case,  it  is  necessary  to  muLtiply  S  s  by  the  joint  probability 


density  function  p(n,U)  and  integrate  with  respect  to  h  ,h  (7) 

X  z 

facilitate  the  evaluation  of  (A.l),  it  is  rewritten  as  follows, 
2  °° 

PQ 


and  U .  To 


'  Pq.  /  |  1 2 '  i  ,2/  PQ  ,  ,2\  1 

<s  'X7-ixi  ;+.xi  ' ‘  S  >  -r-rrr^  i  y 

1  i  1  n  •  a  ^ 

—oo  y 

•exp  iv  x,  +  iv  zl  dx ,  dz 
*-  x  d  zdjdd 

For  high  frequencies  it  can  be  shown  (Sancer  1969,  Brown  1978),  that 
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(A. 5) 
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Thus  (A. 4)  reduces  to 
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<0  ^>=  ~ 
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>  {x2-!x|2|exPLivxxd+  ivzzdJdxd  dZd 
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j  Ixi 


2L  Q 
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>-i<b!^>'2 
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'|exp[i 


iv  x,+  iv  z.jdx,  dz 
x  d  z  d  d  d 


(A. 6) 


y  y 

For  slightly  rough  surfaces  (perturbation  solution) 


DPll 

|  dpq] 

n‘a 

L  yJ 

n  *a 

-  yJ 
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PQ 


(A. 7) 


n-*-a 


where  I)FC^  is  given  by  (35).  Furthermore,  when  the  stationary  phase  approximations 
are  valid  (see  Section  4) 


KQ1 

n*a 

n‘a 

■  y 

L  yJ 

where 


(A. 8) 


DP^/n-a  I-  is  given  by  (64).  Thus  if  jx)2<<  3  when  either  the  physical 
y|ns 

optics  or  perturbation  approximations  are  valid,  the  second  integral  in  (A. 6)  can  be 
PQ 

neglected  and  <o  >  reduces  to  (48). 

PQ  c 

The  scattering  cross  section  for  the  coherent  fields  <o  >  corresponds  to  the  last 
term  in  (48).  Thus 
2 

•  if  v  (\i  l  < . >  pvn  F 1',  . 

X 


PQ  c  o 
'"<3  X>  =  — ~ 
ttA 


X(v  )<. — __>  exp[iv  x  +iv  z]dxdz 
^  n '  a 

"  y 


k2  A 


dpq 

y  n’a 

y 


sinc(v  L  )sinc(v  L  ) 
XX  z  z 


(A.  9) 


in  which  sinc(u)  =  sin(u)/u  and  the  projection  of  the  rough  surface  on  the  xz 
plane  is  given  by 

A  “ 


2L  •  2L 
X  z 


(A. 10) 
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Thus  for  the  specular  case  v  =  v  = 

X  z 


on  •  ^A 

<oPV  =  -°-2 

s  7r 


X(v, ,) 


ipQ 


y  n-a 


(A. Ill 


The  total  scattering  cross  section  corresponds  to  the  first  term  in  (A.l).  Thus 


PQ  T  o 
ttA 


PQ 

-S  >X9(v  , — v  )exp[iv  (x-x ' )+  iv  (z-z')Jdxdz  dx'dz' 
z  y  y  x  z 


(A. 12) 


PQ  c 

However,  except  for  very  rough  surfaces,  where  <o  >  ,  (A. 11)  can  be  neglected, 

PQ  T 

the  total  scattering  cross  section  <o  >  is  evaluated  by  summing  the  incoherent 
and  coherent  terras  (A. 6)  and  (A. 11)  respectively  rather  than  by  directly  evaluating 
(A. 12) (Beckmann  and  Spizzichino  1963),  Thus, 

PQ  T  PQ  c  PQ 

<0  «  <a  +  <o  (A. 13) 
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Fig.  1.  Plane  of  incidence,  scattering  plane  and  reference  x.z  plane 
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SCATTERING  CROSS  SECTIONS  FOR  COMPOSITE  RANDOM  SURFACES— 
FULL  WAVE  ANALYSIS 


I 

Abstract 

The  full  wave  approach  to  rough  surface  scattering  is  applied  to  composite 
rough  surfaces.  In  this  work  the  principal  distinguishing  features  of  the 
individual  rough  surface  is  its  correlation  distance.  Thus  this  model  can  be 
applied  to  scattering  by  rough  seas  as  well  as  hilly  terrain. 

It  is  shown  that  the  full  wave  approach  accounts  for  both  specular  scatter  and 
Bragg  scattering.  The  scattering  cross  section  for  the  composite  surface,  with 
two  or  more  roughness  scales,  is  shown  to  be  a  weighted  sum  of  the  scattering 
cross  sections  for  the  individual  rough  surface  heights.  Shadowing  effects  are 
accounted  for  explicitly  in  the  analysis.  The  full  wave  solutions  satisfy 
reciprocity,  duality  and  realizability  relationships  in  electromagnetic  theory. 


mod  e 1 s 


1.  Introduction 
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Solutions  have  been  derived  for  the  scattering  cross  sections  per  unit  area 

of  the  rough  surface  using  the  full  wave  approach  (Bahar  1981a, b).  For  slightly 

rough  surfaces  the  full  wave  solutions  are  shown  to  reduce  to  the  perturbation 

solution  (Barrick  and  Peake  1968,  Barrick  1970).  When  the  major  contributions 

to  tue  scattered  fields  come  from  the  neighborhood  of  stationary  phase  (specular,) 

points  of  the  rough  surface,  the  full  wave  solutions  are  shown  to  reduce  to  physical 

optics — solutions  (Beckmann  and  Spizzichino,  1963).  Since  the  results  of  the  two 

general  approaches  to  random  rough  scat  ter ing,  per turbat ion  and  physical  optics, 

are  shown  to  be  special  cases  of  the  full  wave  solutions,  the  limitations  of  these 

approaches  can  be  precisely  determined  and  the  differences  between  them  reconciled, 
(Bahar  1981a, b). 

The  principal  motivations  for  this  work  arc  to  extend  the  full  wave  analysis  to 

composite  rough  surfaces  with  multiple  roughness  scales  and  to  explicitly  account 

for  shadowing  effects  in  the  results.  The  main  distinguishing  feature  of  the 

individual  rough  surface  h.  is  its  correlation  distance  i  .  However,  no  res- 

i  i 

trictions  are  made  on  the  variance  of  the  rough  surface  heights  o  ,.  This  work  can 

oi 

therefore  be  applied  to  scattering  by  rough  seas  or  by  hilly  terrain.  In  tht  treatment 
of  composite  rough  surfaces  by  Brown  (1978)  the  feature  that  distinguishes  the  two 
surfaces  considered and  h^is  the  surface  wave  number  k^  where  the  "spectral 
splitting"  occurs. 

The  principal  expressions  derived  for  the  scattering  cross  section,  using 
the  full  wave  approach,  are  summarized  in  Section  2.  The  high  frequency  approxima¬ 
tions  arc  extended  in  Section  3,  to  cases  in  which  the  reference  planes  of 
incidence  and  scatter  are  not  coplanar.  In  Section  4  a  composite  model  of  the 
rough  surface  with  different  roughness  scales  Is  analyzed.  A  two  scale  model  is 
first  considered  and  the  result  Is  given  by  equation  (44).  Using  this  model, 
the  full  wave  analysis  is  shown  to  account  for  both  specular  scatter  and  Bragg 
scattering  (54).  .’oi  composite  surfaces  with  N  uncorr elated  surface  heights 
h^,  the  solution  is  given  by  (55).  In  general, it  is  shown  that  the  scattering 
cross  section  for  the  composite  rough  surface  is  a  weighted  sum  of  the  cross  sections 
for  the  individual  rough  surfaces  h^.  In  Section  5  shadowing  effects  are  explicitly 


2.1  Formulation  of  the  Problem 


The  expression  for  the  scattering  cross  section  per  unit  area,  based  on  the 


full  wave  solutions  for  the  incoherent  fields,  is  given  by  (Bahar  1981b) 
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y  y 


jexpUvxxd  +  ivzzdJdxd  dzd 
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The  superscripts  P,Q  =  V,H  denote  vertically  and  horizontally  polarized  waves  respec¬ 
tively.  The  first  superscript  denotes  the  polarization  of  the  scattered  wave  and  the 
second  the  polarization  of  the  incident  wave.  The  symbol  <  >  denotes  statistical  average 
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in  which  D  x(r)  the  scattering  coef f icients  ,U (r)  is  the  shadow  function 

and  n  is  the  unit  vector  normal  to  the  rough  surface  (Bahar  1981b).  The  vector 

v  is 

v  =  k  (n^-n*)  «*  v  a  +v  a  +  v  a  (3) 

o  xyyyzz  v 

The  free  space  wave  number  is  k^  and  the  unit  vectors  in  the  direction  of  the 

incident  and  scattered  waves  are  n^  and  respectively.  The  rough  surface  characteristic 

function  and  the  joint  characteristic  function  are  x(Vy)  ar>d  X2^Vy,_Vy^  respectively. 

It  is  assumed  that  rough  surfaces  are  statistically  homogeneous  and  isotropic,  thus 

the  surface  height  correlation  function  <  h (x ,z)h (x ' , z ' )  >  depends  only  on  the 

distances  (see  Fig.  1) 

|  rdj=  (x-x’)ax  +  (z-z')az  =  xd  ax  +  7d  az  ,  (4) 

and  the  correlation  distance  is  smaller  than  the  width  of  the  illuminated  surface. 
Furthermore  in  (1)  it  is  assumed  that  the  rough  surface  height  and  slope  are 
uncorrelated.  It  is  also  assumed  that  distribution  of  the  shadow  function  lUr^) 

depends  primarily  opon  the  probability  density  function  for  the  slope  of  the 
rough  surface  (Sancer  1969).  The  elements  DP^(P,Q=V,H)  of  the  2x2  matrix  D  are 


given  by  (Bahar  1981a, b) 
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in  which,  for  the  incident  and  scattered  waves  (denoted  by  superscripts  i  and  i 
respectively) 
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The  transformation  matrices  and  T  relate  the  vertically  and  horizontally 
polarized  components  with  respect  to  the  reference  plane  (normal  to  a^)  to  the 
vertically  and  horizontally  polarized  components  with  respect  to  the  local  plane 
normal  to  the  unit  vector 


n  =  n(  ,h)=Vf/|Vf]  =  (-h  a  +  a  -h  a) / [h4  +  1  +  h^ | 
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in  which 


The  function 


h  =  3h/c)x  and  h  =  9h/9z 
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f  (x,y,z)  =  y  -  h(x,z  ) 


(7b) 
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defines  the  rough  surface  between  medium  0,  y  >  h(x,z),  and  medium  1 ,  y  <  h(x,z).  Th 
angle  between  the  reference  and  local  planes  of  incidence  and  scattering  are  p1  and 
respectively.  The  elements  of  the 
matrix 
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are  functions  of  the  angles  of  incidence  and  scatter  in  the  local  coordinate  system 
and  the  relative  permittivity  c  and  permeability  of  the  two  media  surrounding 

the  interface  f(x,y,z)  -  0.  Thus  FPQ(i’«V,H)  are  explicitly  functions  of  the  unit 
vectors  (see  Bahar  1981a, b) 
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and  the  electromagnetic  parameters  and  respectively  (Bahar  1981b). 

It  has  been  shown  (Bahar  1981b),  that  on  assuming 

k2  <  h2  >  =  k2  a2  «  1,  n*a  -  1  and  U(r  )  -  1  (12; 

o  o  o  y  s 

the  full  wave  solutions  for  the  scattering  cross  sections  per  unit  area  reduce  to 
the  first  order  perturbation  solutions  (Barrick  and  Peake  1968,  Barrick  1970) 
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in  which 
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The  upper  and  lower  signs  in  (L4a)  are  used  for  Q=H  and  V  respectively.  The 
spectral  density  of  the  rough  surface  height  function  is  given  by  the  Fourier 
transform  of  the  height  correlation  function  (Barrick  1970,  Ishimaru  1978). 
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Subject  to  the  perturbation  approximation,  the  local  angles  of  incidence  and  scatter 

are  approximated  by  the  angles  of  incidence  and  scatter  with  respect  to  the 

reference  plane  normal  to  a^.  Thus  for  instance 

„in  .in  -i  -  -i  -  ,i  _i  /-.c-, 

C  “  cosO  =  -n  *n  2  -n  *a  =  cosd  =  C  .  (15) 
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Under  special  conditions  (Bahar  1980b),  the  high  frequency  approximation  of  the 

full  wave  solution  (1)  is  obtained  by  setting 

n  =  n  ,  | n1  a  n]=[n^a  n  J  =  In*  a  n^,  =  0  and  U(r  )  =  1  (16) 
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Thus,  if  the  planes  of  incidence  and  scatter  with  respect  to  the  reference  plane 
(y^o)  and  the  local  plane  at  the  specular  points  (where  n  »  n  )  are  coplanar,  and 
if  the  specular  points  exist  on  the  rough  surface. 
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in  which  -ff  is  given  by 
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3nd  R,  (0  )  is  the  Fresnel  reflection  coefficient  for  vertically  and  horizontally 
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polarized  waves  (P=V,H)  evaluated  at  stationary  phase  points  where 

n  *■  n  =  (n^-n^)/jn^  -n'M 
Thus  S 
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(n**n)^  =  (n^*n)“  =  cos^Os  ^  (CS)^=1^[l+cos0^cos0^-sin4)^sin4>^cos((f^-v^)  j  (  19) 
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With  the  exception  of  the  reflection  coefficient  R,  (0’’),  which  is  evaluated 

lo  o 

at  the  stationary  phase  (specular)  point  rather  than  at  the  angle  of  incidence 

9^,  the  solution  (17  )  is  the  same  as  the  Physical  Optics  solution(£arrick  1970, beck.  . 

anp  Srizzicnino  1963,  Ishimaru  1978).  Thus  for  very  good  conducting  boundaries 
P  s  PP 

(|R  (0  )|  1),  <  Ow  >  is  independent  of  polarization  in  the  Physical  Optics 

limit.  Furthermore  in  view  of  the  Kronecker  delta  <5^  in  (17)  there  is  no 
PQ 

depolarization  (<  >  =  0  for  Pj^Q) .  This  is  due  to  the  assumption 

Ln*  n^]  *  0  (  16) . 

Using  the  full  wave  solution  (1)  as  a  starting  point  for  the  present  analysis, 
the  scattering  cross  sections  are  derived  for  cases  in  which  neither  the 
assumptions  made  for  the  perturbation  analysis  (  12)  nor  the  assumptions  made 
for  the  Physical  Optics  analysis  (  16)  are  valid.  For  instance,  even  for  very 
high  frequencies,  the  assumption  that  the  major  contribution  to  the  (incoherent) 
scattering  cross  sect ’on  come  from  the  stationary  phase  (specular)  points  of 
the  rough  surface  cannot  in  general  be  satisfied  (Bahar  198  lb). 

In  view  of  experimental  evidence  that  neither  the  perturbation  nor  the 
Physical  Optics  approach  aptly  determine  the  scattering  cross  section  for  all  ii*  and 
n  ,and  the  realization  that  simple,  single  scale  models  of  rough  surfaces  are  not 
suitable  for  a  large  variety  of  relevant  problems,  expressions  for  the  scattering  cross 
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sections  have  been  derived  for  composite  rough  surfaces  (Barrick  and  Peake  1968, 
Wright  1968,  Valenzuela  1978,  Brown  1978,  Burrows  1973).  The  full  wave  approach 
will  therefore  also  be  applied  to  composite  rough  surfaces  and  compared  with  earlier 
approaches  to  rhis  problem. 

3.  High  Frequency  Approximation  When  the  Reference 
Planes  of  Incidence  and  Scatter  Are  Not  Coplanar. 

The  Physical  Optics  solution  (17)  was  derived  on  assuming  tha-  the  principal 
contributions  to  the  incoherent  scattered  field  come  from  the  neighborhood  of 
stationary  phase  (specular)  points  of  the  rough  surface  where  n  =  n  .  Furthermore 
it  is  assumed  that  the  phase  in  the  integrand  of  the  full  wave  solutions  varies 
very  rapidly  fcr  the  shadow  region  and  therefore  U(r  )  can  be  set  equal  to  unity 
(Bahar  1981 b).  These  assumptions  are  satisfied  if  the  angle  ~ -9*  is  larger  than 
the  mean  value  of  the  slope  of  the  rough  surface  8q  (Bahar  198lb)  .  Thus  to  satisfy 
these  assumptions  it  is  necessary  that 

tan^f-O1’^  cot0i’^>tan  8  =  2o/£.  (20) 

Z  o  o  o 

where  o  and  l  are  the  variance  and  the  correlation  distance  for  the  rough 
surface.  However,  the  additional  assumption  that  the  reference  planes  of 
incidence  and  scatter  are  coplanar 

[,ni  a^  n^]  =  0  (21) 

is  not  satisfied  in  general  even  when  (20)  is  satisfied.  While  the  assumption 
(21)  simplifies  the  Physical  Optics  solution  (17)  considerably,  it  results 
in  no  depolarization.  In  this  section  the  restriction  (21)  is  lifted  and  the 
appropriate  high  frequency  results  are  derived. 

The  value  of  the  matrix  D  at  the  stationary  phase  points  is  obtained  by 


setting  n  "  ng  (18)  in  equation  (5).  Thus 
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The  angle  93  is  given  by  (  19)  and 
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Thus  the  stationary  phase,  high  frequency  approximation,  for  the  scatterir 
cross  sections  for  P  »  V,H  and  Q  =»  11, V  (P/Q)  are 
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For  highly  conducting  boundaries 
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For  backscatter  however,  since 
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both  (17)  and  (26)  reduce  to 
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For  slightly  rough  surfaces  (Bahar  1981b) 
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Thus,  except  for  normal  incidence  the  perturbation  and  stationary  phase  approximations 
PO 


ror  <o 


for  slightly  rough  surfaces  are  not  in  agreement.  However  it  should  be  noted 


that  for  backscatter  both  the  small  slope  perturbation  solution  n  =  a  and  stationary 

PQ  ^  PQ 

phase  appeoxiraat ions  for  <o  S  (P*Q)  vanish.  The  full  wave  solution  for  <o  >g  (P*Q) 

PO  ” 

does  not  vanish  since  S  ‘  (1)  vanishes  for  backscatter  only  at  the  specular  points  or 
where  fi  ~  ay . 

A.  Composite  Surface  with  Multiple  Roughness  Scales 

In  this  section  it  is  assumed  that  the  composite  surface  under  consideration 

has  two  or  more  distinct  classes  of  roughness  that  are  uncorrelated.  Consider 

first  the  case  in  which  Che  two  statistically  independent  surface  height  variances 

and  a  „  are  small  and 
o  2 
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Thus  i  can  be  factored  out  of  the  integral  Q)  and  since 
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in  which  the  surface  height  correiation  functions  <  h^h^  >  and  <  h^h^  -•  are 
expressed  in  terms  of  the  correiation  coefficients  and  that  depend  on); 
upon  distance  r^  (4).  Thus  the  integral^'  (17b)  reduces  to 
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The  characteristic  functions  and  the  spectral  density  functions  for  the  rough 
surface  height  h^  are  and  respectively.  Thus, in  thi"  case  since  in 

Q)  1  n  ~  a^,the  scattering  cross  sections  for  the  incoherent  fields  is  the  sum 
of  the  scattering  cross  sections  for  the  individual  slightly  rough  surfaces 
and  h^. 

For  the  more  general  case  the  surface  height  variances  o  ,  and  a  are  not 

ol  o2 

assumed  to  be  smalJ .  However,  the  correlation  distance  is  assumed  to  bo  much 

smaller  than  Thus  by  definition 


Cf  (p=£^)  =  C^/e  and  *=  C^/e 


in  which  p  «*  (x2+z2)1;i,e  is  the  Neperian  number  and  it  is  assumed  that, 
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where  9.^  is  a  constant,  in  this  case  <U  >  (  1)  is  given  by 
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The  significance  of  the  assumption  (37)  is  that  for  distances  0  <  p  <  £j2 

X^O)  =  1  >  X2  >  I X1 1 2  x\  *  -  1 


(39) 
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and  for  distance  p 
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Assuming  for  example  normal  distributions  for  h. 
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and  the  scattering  cross  section  can  be  written  as 
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and  .SP  >.  the  statistical  average  of  (2)  is  given  by 

f  SPQ  p (n,U)dh  dh  dU 
J  x  z 

where  p(n,U)  is  the  joint  probability  density  function  (Sancer  1969).  If  the  variance 
of  the  surface  height  h^  is  small^lx^j  -  l,and  (43)  reduces  to 
PQ  PO  ,  PQ  2  2 

■•o  >  =  <0  >,  +  <o  ■>  ,  v  o  .  <<  1  (43) 

1  2  y  ol 

PO 

Before  the  transform  in  (44b)  can  be  evaluated  it  is  necessary  to  determine  <S  > , 
(44c)  for  the  composite  rough  surface.  When  the  stationary  phase  approximations 
are  valid  for  high  frequencies  and  shadowing  effects  are  not  very  significant  n 
can  be  replaced  by  n^  in  (45)  and  as  jn  (26) 

'  -NPQ 


PQ 


(42) 


(43) 
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When  the  mean  value  of  the  slope  of  the  composite  surface,  Bq  (20)  is  very  small 


tanB 
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(47) 


rr-ii 


the  perturbation  approximation  for  (44c)  can  be  used  as  in  (13) 

<sPQ>  *  sP(^  =  |dp<5|2 


For  near  grazing  angles  where 

cotOi,E  «  can  ft  =  la'l  (49) 

o  o 

shadowing  effects  become  significant  and  even  at  high  frequencies  (46)  cannot  be 
used.  In  this  case,since  U (r  )  =■  1  only  for  those  portions  of  the  rough  surface 


that  are  almost  horizontal  (n*a^  =  1),  the 


perturbation  approx ima t ion 


n  ay  could  be  used  for  near  grazing  angles  (49)  even  at  high  frequencies. 

For  the  general  case  however,  it  is  necessary  to  evaluate  <SP(^>  (44c),  using  the 

statistics  of  the  slopes  of  the  composite  surface  (see  Sec.  5).  Provided  that  (37) 

satisfied  and  either  the  perturbation  or  stationary  phase  approximations  are 
PO 

valid  such  that S  x  is  independent  of  position,  (46)  ,  (48)  ,  the  scattering 

cross  sections  (1)  can  be  expressed  as  follows: 
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For  rough  surfaces  with  normal  distributions  (Beckmann  and  Spizzichino  1963) 

Tr(v  J  . £. )2exp(-v2u2  -v2  S.?/4)  =  tt2v2  W(v  ,v  ) ,  (v  a  ,)2  <<  1 
y  oi  i  r  y  oi  xz  i  y  x  z.  y  oi 

oo  (v  O  )  2R1 

4} .  =  {nf.2exp(-v2a2  - - exp(-v2  f,2/4),  v  a  .  =  1 

\  y  oi  *•  m'm  r  xz  i  y  oi 


Uf.  r  /  „  .  ,2 

i  ,■?'  expL-(v  v  a  . ) 

2  2  r  xz  y  oi 

tf  rr  J 


(v  n  )  »  i 

y  oi 


where 


(vx  +  Vz^ 


The  expression  (50)  simplifies  considerably  if  the  variance  of  the  surface  height 

2  2 

h,  Is  small  (44)  while  the  variance  of  the  surface  height  h,,  is  large  (v  <<  1 
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and  >>  1).  In  this  special  case  the  contribution  f  r  ,n  the  term  <n'L  ^  ■  is 

v  o,'.  2 


dominant  when 


tan  6  =  2o/? 

o 


On  the  other  hand  only  horizontal  portions  of  the  rough  surface  (n  -  a  )  are  both 

y 

illuminated  and  visible  for  grazing  angles,  (49)  (  and  the  contribution  from  tiie 
PO 

term  •  '>  is  dominant  near  grazing  angles.  Thus  if  in  addition  to  condition  (37) 

2  2 

v  a  „  ■  1  and  v  a  '>  1,  at  high  frequencies  (c.n\  can  be  approximated  by 

y  o2  y  v  1 


.  PQ.  o  „PQ  v  PQ  ,  .  _  PQ  PQ 

•  d  v  -  -  —  Ls  x  f.  +  s  y  .jL  '  =  <a  v>,  +  <a 

7.  o  >  1  <*>  .J2‘  1 


PQ  PQ 

in  which  Sa  and  are  given  by  (46)  and  (48)  respectively.  The  scattering 
cross  sections  for  such  composite  surfaces  are  approximately  equal  to  the  sum  of 

PQ  PQ 

the  individual  cross  sections  <0  and  <0  ^>2  for  the  rough  surface  heights 

and  h^  respectively  (Barrick  and  Peake  1968,  Barrick  1970). 

The  result  (43)  can  be  generalized  for  a  composite  surface  represented  by 

the  superposition  of  N  uncorrelated  surface  heights  h  (n  =  1.2...N)  for  which 

n 

the  correlation  distances  SL  satisfy 

n  J 

ll  'K  *12  <<  *2  <<  *23  <<  *3  ••‘•Vi  ^  *N-1,N  <<  *N  (55 

In  this  case 


PQ  y  tt  1  iu-1  1 2  PQ  y  PQ 

<(7  '  *  l  II  IX  I  <0  =  l  W  so 


n  =  l 


n  *- ,  n  n 

n=l 


where 


- §■  |  <SP(W  “lx"!  Jexp[  iv  x  +  iv  z  jdx  dz 


z  ~dJ  d  d 


The  characteristic  functions  for  the  surface  height  h  are  Y»(v  ,-v  )  and 

n  ^2  y  y 

X  (v  ) ,  and  x°  -  !•  In  (53)  the  contribution  <oP^>n  t0  the  total  scattering 
PQ 

cross  section,  >,  due  to  the  rough  surface  h^  is  weighted  by  the  product 


vn  -  fl'|xm~1|2 

m-1 


Since  |x*|  <  1  •  the  weighting  factor  w  is  in  general  less  than  unity. 


(57) 
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It  decreases  monotonically  as  n  increases  even  when  the  variance  o  of 

ox 

2  2 

individual  rough  surfaces  are  small  (v  a  «  1,  i  =  1.2...N). 

y  oi 

The  effects  of  shadowing  have  not  been  accounted  for  explicitly  in 
solution  (54).  To  do  so  it  is  necessary  to  determine  the  statistical 

PQ 

<S  >  (45)  using  the  joint  probabilty  density  function  p(n,U)  for  the 

surface  slopes  Vf  (7)  and  the  shadowing  function  U(r  ).  Both  terms  of 
nroduc  t 


t  he 

the 

averages 

rough 

the 


<SPQ>  [Xo  *  I Xn 1 2 ]  =  K(xd,zd)  ( 58) 

are  functions  of  distance  r.  (4).  The  two  dimensional  Fourier  transform  of  the 

d 

product  K(xd,zd),  (56)  can  be  evaluated  by  determining  the  two  dimensional 

convolutions  of  the  two  dimensional  Fourier  transforms  for  <SP^>  and  ["  v”  “lx” !  ""  .  (Brown 

1978,  1980).  However  for  high  frequencies  the  integral  expression  for  the 

scattered  fields  (Bahar  1981b)  may  also  be  evaluated  at  the  stationary  phase 

points  before  squaring  and  averaging  (Kodis  1966).  Barrick  (1968)  shows 

that  the  results  are  the  same  regardless  of  the  order  of  evaluating  the  integral 

and  the  statistical  average.  This  latter  approach  which  simplifies  the 

analysis  will  be  used  in  Section  5  to  explicitly  account  for  shadowing  in  the 
PO 

evaluation  of  <0  x>  (56). 

n 

5.  Scattering  Cross  Sections  for  Composite  Rough  Surfaces 

When  Shadowing  Effects  are  Accounted  for  Explicitly 
PO 

To  evaluate  <S  (44c)  ■  for  high  frequencies  it  is  convenient  to  represent 
the  joint  probability  density  function  p(n,U)  in  terms  of  the  conditional 
density  P(Ujn)  (Sancer  1969,  Brown  1978,  1980) 


p(n,U)  =  p(n)p(u|n)  (59) 

in  which  p(n)  is  the  density  function  of  the  gradient  of  the  rough  surface 
Vf  -  njVf]  (7)  and, 

p(u|n)  =  ,n*  ]  n)<5  (U-l)  +  !  l-P^  (n^  ,n^  |  n)  ]5  (U) 


(60) 
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The  Dirac  delta  function  is  6  (U )  and  P^Cn^.n^n)  is  the  probability  that  a  point 

will  be  both  illuminated  by  a  source  at  -r  n1  and  visible  at  the  observation  point 

o 

rnf,  given  the  value  of  the  gradient  ii(hx,hz).  Thus  for  high  frequencies  <S 
can  be  replaced  by  its  stationary  phase  approximation  (Sancer  1969) 

,2 


<sPQ> 


<sp3>  = 


D 


PQ 


p(n)p(u|ns)U2(?s)dhx  dhz  dU 


?2 (n^ , ni | n  ) 


(61) 


y'n 


The  expression  for  the  prooability  P_(n^,n*|n  )has  been  given  by  Sancer  (1969) 

£  S 

for  rough  surfaces  with  normal  distributions.  The  above  approximation  for 
PO 

<S  >  is  appropriate  for  scattering  from  the  rough  surface  h2  with  th<  large 
correlation  distance  Thus 


PQ  o 

<a  >,  =  — 

l  TT 


n*a 


P2(nf  2 


(62) 


in  which 


DPQ/n-a 


is  given  by  (22)  and  -f  ^  (35b)  is  given  by  (51)  for 


surfaces  with  normal  distributions.  Thus  for  surfaces  with  (v^o^)  >:>  ^ 2 
proportional  to  the  joint  probability  density  function  P(ng)  for  the  gradient  at 
the  stationary  phase  points  (51). 

For  angles  0i’^  smaller  than  the  mean  value  of  the  slope  8  (53)  , 

o  o 

PQ 

<u  x>2  is  the  dominant  term  in  the  solution  (44).  However  for  angles  larger 
than  B0  (near  grazing  angles  0*’^)  (49),  scattering  due  to  the  small  scale  rough¬ 

ness  dominates.  This  is  either  because  the  joint  probability  density  function  at  the 
stationary  phase  points,  p (ng)  is  very  small  or  because  the  probability  is  very 


small  that  these  stationary  phase  points  are  illuminated  and  visible, 


1 

1 
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( Pjfu*  ,ni  |  n  y  ■  *'  1).  S<w- wring  duo  to  the  smalt  scale*  ^jghness  is  roost 

signi Meant  lor  near  grazing  angles  whore  the  contributions  from  the  specular 

points  are  very  small.  Thus  the  contribution  to  the  total  scattering  cross  section 
PO 

due  to  the  term  -n  (44)  is  primarily  from  the  near  horizontal  portions  of 

the  rough  surface  that  are  both  illuminated  and  visible.  Hence  for  the  term 
PO  PQ 

' o  <S  >  can  be  approximated  by 

<sp«>  - 

O 

y 


|  ppQ 

I  n*a 

)  i  v 


p(n)p(U|n)U  (r  )dhx  dh^  dU 


P  (nf,nl,U-l) 


P.  (n  ,n  tli»l)  °  |P  (n  ,n  jn)p(n)dh  dh 
l  j  Z  *  y 

is  the  probability  that  the  rough  surface  is  both  illuminated  and  visible. 

Thus  for  the  surface  with  the  small  roughness  scale 


n  kV|DP«l 


Pj.  (n*  *nA  ,U«=1)W1  (vx,v2) 

n-a 

y 


in  which  W^(v^,v  )  is  the  spectral  density  function  for  the  rough  surface 


The  scattering  cross  section  for  the  composite  rough  surface  is  therefore 

obtained  by  summing  (65)  and  (62)  while  accounting  for  the  weighting  function 
x  2 

w^  *>  | y  |  in  (56).  This  weighting  function  (w^  <  1)  has  the  effect  of 
damping  the  >'ontr ibuti on  from  the  specuLar  poll  *  because  of  the  superposition  of 
the  surface  (with  the  small  roughness  scale)  on  the  surface  (with  the 
large  roughness  scale). 

On  deriving  the  expression  for  the  scattering  cross  section  an  integral 


containing  the  term 


,  I  h>  vu I  1  D‘wu  . 

A  -  <  <  r-r-  P  -  \  1* 

\  ln'ayi  I  n'ay  j  I  I 


has  been  neglected  (Rahar  1981b).  This  term  can  be  Ignored  either  If  |xl  <<  l 
or  when  the  perturbation  or  stationary  phase  approximations  are  valid. 
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6.  Concluding  Remarks 

Using  the  full  wave  approach,  a  general  expression  has  been  derived  for  Che 
scaccering  cross  seccions  for  composiCe  rough  surfaces  comprised  of  a  superposic ion 
of  N  uncorrelaced  rough  surface  heighcs  Ik.  The  distinguishing  feaCure  of  the  differ¬ 
ent  rough  surface  heights  is  Che  correlation  distance  (55)  and  not  the  surface 
height  variance.  It  is  shown  that  the  scattering  cross  sections  for  the  composite 
surface  are  a  weighted  sum  of  the  scattering  cross  sections  for  the  individual  rough 
surfaces.  Thus  for  a  composite  surface  with  two  roughness  scales,  the  first  sLightly 
rough  and  the  second  very  rough  and  with  the  longer  correlation  distance,  the  scattering 
cross  section  accounts  for  both  Bragg  scattering  and  specular  scatter.  However,  the 
contribution  due  to  specular  scatter,  by  the  very  rough  surface,  is  slightly  damped 
as  a  result  of  the  superimposed  slightly  rough  surface.  Shadowing  effects  have  been 
accounted  for  explicitly  in  the  analysis.  If  the  correlation  distance  JL  is  noL  the 
distinguishing  feature  of  the  different  rough  surface  heights  Ik,  the  starting  point 

of  the  analysis  for  the  composite  surface  is  (38).  Since  (43)  assumes  that  the 
correlation  distance  is  the  distinguishing  feature  of  the  different  rough  surface 
heights  h.  (37),  it  cannot  be  used  in  general, and  the  cross  section  for  the  composite 
surface  is  not  given  by  a  superposition  of  the  cross  sections  for  the  surface 
heights  h^  (44),  (45)  or  (56). 
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